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Abstract 
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^ ' Explicit expressions for the non-singlet and singlet structure functions g\ at 

q ■ the small x-region are obtained. They include the total resummation of double- 

logarithmic contributions and accounting for the running QCD coupling effects. We 
predict that asymptotically the singlet g\ ~ x _As , with the intercept As = 0.86, 
which is approximately twice larger than the non-singlet intercept Ajvs = 0.4. The 
> i impact of the initial quark and gluon densities on the sign of* at x « 1 is discussed. 
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' 1 Introduction 

p , 

Deep inelastic scattering (DIS) is one of the basic processes for investigating the structure 
of hadrons. As is well known, all information about the hadrons participating into DIS 
comes from the hadronic tensor Wp V . The imaginary part of Wp V is proportional to the 
forward Compton amplitude when the deeply off-shell photon with virtuality q 2 scatters 
off an on-shell hadron with momentum p. For the electron-hadron DIS, the spin-dependent 



part,W^, of is 
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S P 9i + (S p p p )g 2 ~ le^txp Sh 1 + S^fa + g 2 ) 



pq 



pq 



where m is the hadron mass, S'jJ and Sp are the longitudinal and transverse (with re- 
spect to the plane formed by p and q) components of the hadron spin S p ; g x and g 2 are 
the spin structure functions. Both of them depend on x = —q 2 /2pq, < x < 1 and 
Q 2 = —q 2 > 0. Obviously, small x corresponds to s = (p + q) 2 ~ 2pq ^> Q 2 . In this 
case, S p \ ~ Pp/m and therefore the part of W s related to g\ does not depend on m. Then 
if Q 2 3> m 2 , one can assume the factorization and regard W* as a convolution of two 
objects (see Fig. 1). The first of them is the probability $ ($ = $ q in Fig. la and $ = $ 3 
in Fig. lb) to find a polarized parton (a quark or a gluon) in the hadron. The second one 
is the partonic tensor W* defined and parametrized similarly to W*. Whereas $ 9iS are 
essentially non-perturbative objects, the partonic tensor W^ u , i.e. the partonic structure 
functions g\ and g 2 , can be studied within perturbative QCD. THe lack of knowledge of 
$ is usually compensated by introducing initial parton distributions that are found from 
pheneomenological considerations. On the contrary, there are regular perturbative meth- 
ods for calculating the structure functions in the partonic tensor W*. The best known 
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Figure 1: Representation of the hadronic tensor as the covolution of the fragmentation 
functions and the partonic tensor. 

instrument to calculate the structure functions to all orders in a s is the DGLAP[1] ap- 
proach. Once applied to the description of the experimental data, DGLAP provides good 
result s[2]. The extrapolation into the small-x region of DGLAP predicts an asymptoti- 
cal behaviour ~ exp(-y/C ln(l/x) In In Q 2 ) for all DIS structure functions (with different 
factors C). However, from a theoretical point of view, such an extrapolation is rather 
doubtful. In particular, it neglects in a systemetical way contributions ~ {a s \v?{l/ x)) k 
which are small when x ~ 1 but become large when x <C 1. The total resummation 
of these double-logarithmic (DL) contributions made in Refs. [3] and Ref. [4] for the 
non-singlet (gi S ) and singlet g\ respectively leads to the Regge (power-like) asymptotics 
9i{9\ S ) ~ (1/ X ) A ((^/^) Ans ), with A DL ,A®g being the intercepts calculated ib the 
double-logarithmic approximation (DLA). The weakest point of Refs. [3], [4] is the fact 
that they keep a s fixed (at some unknown scale). It leads therefore to the intercepts 
A DL , A^s explicitly depending on this unknown coupling, whereas a s is well-known to 
be running. The results of Refs. [3], [4] led many authors (see e.g. [5]) to suggest that the 
DGLAP parametrization a s = a s (Q 2 ) has to be used. However, according to results of 
Ref. [6], such a parametrization is correct only for x ~ 1 and cannot be used for x <C 1. 
The explicit dependence of a s suggested in Ref. [6] has been used to calculate both g^ s 
and g\ at small x in Refs. [7]. The present talk is based on the results obtained in those 
papers. 

Instead of a direct study of g±, it is more convenient to consider the forward Compton 
amplitude A related to gi as follows: 

gi (x,Q 2 ) = -^ s A(s,Q 2 ) . (2) 

7T 

We cannot use DGLAP for studying g l or A at small x because it does not account 
for double-log (and single- log) contributions which are independent of Q 2 . In order to 
account for the double- logs of both x and Q 2 , we need to construct two-dimensional 
evolution equations that would combine the x- and Q 2 - evolutions. On the other hand, 
these equations should sum up the contributions of the Feynman graphs involved to all 
orders in a s . Some of those graphs have either ultraviolet or infrared (IR) divergencies. 
The ultraviolet divergencies are regulated by the usual renormalization procedure. In 
order to regulate the IR ones, we have to introduce an IR cut-off. We use the IR cut-off 
fi in the transverse momentum space for momenta ki of all virtual quarks and gluons: 



2 



fi < k i± 



(3) 



where ki± stands for the transverse (with respect to the plane formed by the external 
momenta p and q) component of fcj. This way of regulating the IR divergencies was 
suggested by L.N. Lipatov and used first in Ref. [8] for quark-quark scattering. Using this 
cut-off p, A acquires a dependence on p. Therefore, one can evolve A with respect to /i, 
constructing thereby some Infrared Evolution Equations (IREE). As A = A{s/ p 2 , Q 2 /p 2 ), 

-p 2 dA/dp 2 = OA/ dp + dA/dy (4) 

where p = ln(s/p 2 ) and y = ln(Q 2 /p 2 ). Eq. (4) is the lhs of the IREE for A. In order to 
write the rhs of the IREE, it is convenient to use the Sommerfeld- Watson transform 

^(s/p 2 n(u;)F(u,Q 2 ) (5) 
-too Ztti 

where £(cu) is the negative signature factor, £(c<j) = [1 — e~ m0J ]/2 ittuj/2. It must be 
noted that the transform inverse to Eq. (5) involves the imaginary parts of A: 

2 r°° 

F(u, Q 2 ) = — dpe~^A(s, Q 2 ) . (6) 

1TLU JO 

Notice that contrary to the amplitude A, the structure function g\ does not have any 
signature and therefore t;(u) = 1 when the transform (5) is applied directly to g±. 



2 Infrared evolution equations for g\ 

When the factorization depicted in Fig. 1 is assumed, the calculation of g\ (we will not use 
the superscript "s" forgx singlet, though we use the notaion g^ s for the non-singlet g\ ) is 
reduced to calculating the Feynman graphs contributing the partonic tensor W^ v depicted 
as the upper blobs in Fig. 1. Both cases, (a), when the virtual photon scatters off the 
nearly on-shell polarized quark, and (b), when the quark is replaced by the polarized gluon, 
should be taken into account. Therefore, in contrast to Eq. (2), we need to introduce two 
Compton amplitudes: A q and A g corresponding to the upper blob in Fig. la and Fig. lb 
respectively. The subscripts "q" and "g" refer to the initial partons. Therefore, 

g 1 (x,Q 2 )=g q (x,Q 2 )+g q (x,Q 2 ), (7) 

where 

g q = -%A q (s, Q 2 ), g g = -%A g (s, Q 2 ) (8) 

TT TT 

Let us now construct the IREE for the amplitudes A qt9 related to g\. To this aim, let 
us consider a virtual parton with minimal k±. We call such a parton the softest one. If it 
is a gluon, its DL contribution can be factorized, i.e. its DL contribution comes from the 
graphs where its propagator is attached to the external lines. As the gluon propagator 
cannot be attached to photons, this case is absent in IREE for A q ^ g . The second option 
is when the softest partons are a t-channel quark-antiquark or gluon pair. It leads us to 
the IREE depicted in Fig. 2. Applying the operator —p 2 d/dp 2 to it, combining the result 
with Eq. (4) and using (5), we arrive at the following system of equations: 
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Figure 2: Infrared evolution equations for the amplitudes A q ,A g . 



d 1 

{co + —)F q {u,y) = — [F qq {u)F q {u,y) + F qg {u)F g {u,y)} , 



(9) 



The amplitudes F q , F g are related to A q , A g through the transform (5). The Mellin ampli- 
tudes Fik, with i,k — q,g, describe the parton-parton forward scattering. They contain 
DL contributions to all orders in a s . We can introduce the new anomalous dimensions 
Hik = (l/87r 2 )Fjfc. The way of writing the subscripts "q,g" corresponds to the DGLAP- 
notations. Solving this system of eqs. and using Eq. (8) leads to 



g q (x, Q 2 
100 duo 



2m 



(1/xY C + (u) 



•°° duj 
too 2m 
X 



(l/xr\C+(uj)e n+y + C-(u)e Q - y 



(10) 



2H, 



</</ 



2H, 
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The unknown factors C±(u) have to be specified. All other factors in Eq. (10) can be 
expressed in therms of H ik : 



a 



X = H gg- H qqi R = ( H gg ~ H qqf + ^ H qg H gq 
1 



Hqq + H gg ±J (Hgg - H, 



L gg) 



^HggHgg 



The anomalous dimension matrix was calculated in Ref. [7]: 



where 



dqq — 



A(u)C F 



2n 



a gg ~ 



hk — dik + Vik, 

2A(lo)N n f A'(w) 



7T 



a gq — 



2ir 



dgg — 



A'(uo)C F 



7Y 



and 



with 



7T 



c r 



m " ~ 2N ) m " 



-2N 2 , m qg =n r 



N 



m r , 



-NC F . 



(13) 
(14) 
(15) 

(16) 
(17) 
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We have used here the notations Cp = 4/3, N = 3 and nf = 4. The factors A and D 
account for running a s . They are given by the following expressions: 



m = I 



v r 

rf + 7T 2 Jo 



dpe 



■UIp 



1 

262 



dpe- up \n((p + ri)/ri) 



(p + ^ + n 2 
P + V 



+ 



(18) 
(19) 



-(p + ^) 2 + 7T 2 p + rji 

with r\ = \n(p 2 / Aq CD ) and b = (33 — 2nj)/127r. A' is defined as A with the % 2 term 
dropped out. Now we can specify the coefficients C± of Eq. (10). When Q 2 = p 2 , 



g q = Aq(x ), g g = Ag(x ) 



(20) 



where Aq(xo) and Ag(xo) are the input distributions of the polarized partons at x = p 2 / s. 
They do not depend on Q 2 . Using Eq. (20) allows to express C±(u) in terms of Aq(u) and 
Ag(u), which are related to Aq(x ) and Ag(x ) through the ordinary Mellin transform. 
Indeed, 



C+ + a 



Aq, 



C + ^R + C-^*=Ag. 



2H, 



19 



2H, 



(21) 



'19 



This leads to the following expressions for g q and g g : 



9q 



(x,Q 2 ) = £2£-(l/xr[(A^Aq + BAg)e n +y + (A^Aq - BAg)e n - y ] , (22) 



g g (x,Q 2 ) 
with 



r° K^0-/ x ) u \( EA Q + A^Ag)e n+y + ( - EAq + A^Ag)e n - y 
J -too 2tti l ^ ' ^ ' 



A^ = ( 1 -±^=), B = 



H, 



'19 



E 



H, 



gq 



(23) 



(24) 



lyjw VR VR 
Eqs. (22, (23) express g\ in terms of the parton distributions Aq{uS) and Ag(u). How- 
ever, they are related to the distributions Aq(x ) and Ag(x ) at very low x: Xq ~ p 2 /s <C 
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1. Therefore, they scarcely can be found from experimental data. It is much more useful 
to express g q , g g in terms of the initial parton densities Sq and Sg defined at x ~ 1. We can 
do it, using the evolution of Aq(x ), Ag(x ) with respect to s. Indeed, the s-evolution 
of Sq, Sq from s ~ yU 2 to s ^> fi 2 at fixed Q 2 (Q 2 = /i 2 ) is equivalent to their x-evolution 
from i ~ 1 to x < 1. In the w-space, the system of IREE for the parton distributions 
looks quite similar to Eqs. (9). However, the eqs for Aq, Ag are algebraic because they 
do not depend on Q 2 : 



Aq(uj) = (< e 2 > /2)Sq(u) + [H qq (u)Aq(u) + H qg (u)Ag(u)} , 

Ag(u) = (<e 2 q > /2)Sg(u) + (l/u)[H gq (oj)Aq(oj) + H gg (oj)Ag(u;)] . (25) 

where < e 2 > is the sum of the quark electric charges (< e 2 >= 10/9 for rif = 4), Sq 
is the sum of the initial quark and antiquark densities and Sg = — (A'(oj)/2ixuu 2 )Sg is the 
starting point of the evolution of the gluon density Sg. It corresponds to Fig. lb where 
the upper blob is substituted by the quark box. Solving Eqs. (25), we obtain: 

Aq = (<e 2 >/2) - H " )5q + uHt »* 9] (26) 

q 1 q 1 } [uj 2 - uj(H qq + H gg ) + (H qq H gg - H qg H gq )} ' 1 > 

Ag = (< e 2 > /2) [^H gq Sq + u(u - H qq )Sg] (2?) 

q [u 2 - uj(H qq + Hgg) + (H qq H gg - H qg H gq )} 

Then Eqs. (22,23,26,27) express g\ in terms of the initial parton densities Sq, Sg. 

When we put H qg = H gq = H gg = and do not sum over e q , we arrive at the expression 
for the non-singlet structure function g^ s : Obviously, in this case = B = E = f2_ = 
0, = 1, f2 + = H qq . However, the nonsinglet anomalous dimension H qq should be 

calculated in the limit b gg = b qg = b gq = 0. We denote such H qq = H NS . The explicit 
expression for it is: 

H NS = (1/2) [uj - ^ 2 - Ab qq ] . (28) 

Therefore, we arrive at 

^-!/_':^(^)(v<(«v,r™. m 

3 Small-x asymptotics for gi 

When x — > and Q 2 ^> ji 2 , one can neglect contributions with in Eqs. (10). As is 
known, g\ ~ (l/x) w ° at x — > 0, with cj being the position of the leading singularity of 
the integrand of g\ . According to Eqs. (12), the leading singularity, tu NS for g^ s is the 
rightmost root of the equation 

cj 2 - Ab qq = (30) 
while the leading singularity, u for g 1 is the rithmost root of 

u 4 - A{b qq + b gg )u 2 + 16(6^6^ - b qg b gq ) = . (31) 

In our approach, all factors b^ depend on rj = ln(/x 2 / Aq C d), so the roots of Eqs. (30,31) 
also depend on rj. This dependence is plotted in Fig. 3 for uj ns and in Fig. 4 for uj . Both 
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Figure 3: Dependence of the intercept u on infrared cut-off r\ = ln(jU 2 / Aq C d) for g[ 
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Figure 4: Dependence on r\ of the rightmost root of Eq. (31), c^o- Curve 2 corresponds 
to the case when gluon contributions only are taken into account; curve 1 is the result of 
accounting for both gluon and quark contributions. 
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the curve in Fig. 3 and the curve 1 in Fig. 4 have a maximum. We denote this maximum 
as the intercept. Therefore, 



/1 \ Ans fo 2 \ Ans/2 /1 \ As fo 2 \ As/2 

(32) 



and we find for the intercepts 

A NS « 0.4, A s « 0.86 (33) 

and Z 1 = —1.2, Z 2 = —0.08. This implies that g^ s is positive when x — > whereas 
</! can be either positive or negative, depending on the relation between 5q and 5g. In 
particular, is positive when 

15 6q + 6g<0. (34) 

otherwise it is negative. In other words, the sign of g\ at small x can be positive if the 
initial gluon density is negative and large. 



4 Conclusion 

The total resummation of the most singular (~ a™/u 2n+1 ) terms in the expressions 
for the anomalous dimensions and the coefficient functions leads to the expressions of 
Eqs. (7,22,23,29) for the singlet and the non-singlet structure functions g\. It garantees 
the Regge (power-like) behaviour (32) of gi, g± when x — > 0, with the intercepts given by 
Eq. (33). The intercepts Ajvs, As are of course obtained with the running QCD coupling 
effects taken into account. The value of the non-singlet intercept Atvs — 0.4 is now con- 
firmed by several independent analyses [10] of experimental data. The value A s = 0.86 of 
the singlet intercept is in a good agreement with the estimate As = 0.88 ± 0.14 obtained 
in Ref. [11] from analysis of the HERMES data. 

Another interesting point to discuss is the sign of these structure functions. Eq. (29) 
states that g^ s is positive both at x ~ 1 and at x <g; I. The situation concerning the 
singlet g>i is more involved: being positive at x ~ 1, the singlet g\ can remain positive at 
a; C 1 only if the initial parton densities obey Eq. (34), otherwise it becomes negative. 
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